RANDOM AND EXHAUSTIVE GENERATION OF PERMUTATIONS 

AND CYCLES 



MARK C. WILSON 

Abstract. In 1986 S. Sattolo introduced a simple algorithm for uniform random gener- 
ation of cyclic permutations on a fixed number of symbols. This algorithm is very similar 
to the standard method for generating a random permutation, but is less well known. 

We consider both methods in a unified way, and discuss their relation with exhaustive 
generation methods. We analyse several random variables associated with the algorithms 
and find their grand probability generating functions, which gives easy access to moments 
and limit laws. 



1. The algorithms 

Basic notation. For each n > 1, we denote by Sn the symmetric group (set of all 
permutations under the operation of composition) on the set [n] := {1, . . . ,n} . Suppose 
that n > 2. For each vr G S„_i, let vr* be its extension to S„, by definition the element 
of S„ that fixes n and agrees with vr on [n — 1]. The map * is injective on each S„. Dually, 
for an element p € Sn that fixes n, is the restriction of p to Sn-i- The map * is 
onto each S„^i and * followed by ^ is the identity on S„_i. Thus we may consider that 
Sn C Sn+i and let S be the union of all S„ (formally, we consider the direct limit induced 
by the natural inclusion maps [n] — > [n + 1]). Each element vr of S belongs to a maximal 
Sn, where n is the largest integer moved by vr; we define n(vr) to be this value of n. 

The action of vr G S on i is denoted by zvr. We use the standard representation as 
words throughout; the element vr of S„ is written as the word Ivr • • • nvr. Let C„ be the 
set of n- cycles of S„ (recall that an element of S is a fc-cycle if and only if its action has 
a single nontrivial orbit, and this orbit has size k). When n = 1, our convention is that 
Cn = S„,. A 2-cycle is called a transposition, and we denote by T(i,j) the transposition 
that exchanges i and j and fixes all other symbols. 

Finally, for vr € S„, we define g(vr) = mr~^, and Q{ir) = rnr = q{7r^^). 

Random generation. Pseudocode for the two algorithms discussed below is shown in 
FigurelH Note that for notational convenience we will consider permutations on {0, . . . , n— 
1} instead of {1, . . . ,n} in that figure, as well as in section [H 

The standard algorithm [Knutl969l 3.4.2, Algorithm P] for uniformly generating a ran- 
dom permutation of [n] is as follows. Start with the identity permutation. There are 
n — 1 steps. At the ith. step, a random position j is chosen uniformly from [n — i + 1] 
and the current element in position j is swapped with the element at position n — i + 1. 
Example: the permutation 25314 G S5 is formed by choosing j = 4,1,3,1 in that order. 
Knuth attributes this algorithm to R. A. Fisher and F. Yates |FY1938] . and a computer 
implementation was given by Durstenfeld [Dursl964] : it is often called the "Fisher- Yates 
shuffle" or t he "Knuth shuffle" . 

S. Sattolo |Sattl986] introduced a very similar algorithm for uniform random generation 
of an element of C„. The only difference in the algorithm is that the possibility j = n — i+1 
is disallowed; j is chosen uniformly at random from [n — i]. 
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algorithm fisher-yates 

Input: positive integer n 
Output: permutation tt G 

begin 

for i from to n — 1 do 

7r[i] := i 
if (n > 1) then 

for k from n — 1 downto 1 do 

j := raiid(O..A;) 

TT <— swap(7r, r, k) 

return(7r) 
end 



algorithm sattolo 

Input: positive integer n 
Output: permutation tt G C„ 

begin 

for i from to n — 1 do 

7r[i] := i 
if (n > 1) then 

for k from n — 1 downto 1 do 

j := raiid(O..A; — 1) 

TT <— swap(7r, r, k) 

return(7r) 
end 



Figure 1: The random generation algorithms 



Algebraic description. In terms of multiplication in the symmetric group S„, the de- 
scription of the Fisher- Yates algorithm is as follows. For each i,j € [n], let T{i,j) denote 
the transposition that exchanges i and j. Then each algorithm starts with the identity 
permutation tt. At each value of k, the statement tt <— 7rr(r7r,_77r) is executed. In terms of 
multiplication on the left, wc have the following. Observe that the transposition r(r7r,j7r) 
equals the conjugation TT^^T{r, j)TT. Thus we can rewrite the step above as tt <— r(r, j)7r. 

Definition 1.1. We call a product ti • • • of transpositions a triangular product if 
for each i < k, Ti exchanges i + 1 and some j < i + 1. The product is called a strict 
triangular product if always j <i. 

It follows directly from the above discussion that each execution of the Fisher- Yates (re- 
spectively Sattolo's) algorithm yields a triangular (respectively strictly triangular) product 
of ri — 1 transpositions in Sji- Furthermore these maps are 1—1. To see this, note that given 
TT we can reconstruct for each i. This is because fixes all elements greater than i + 
so that n • r„_i = n • tt. This determines the transposition r„_i, and the result follows by 
induction on n. 

We may therefore define maps ^ : 7r„_i x [ra — 1] ^ 7r„ and {i,q) : 7r„ — >■ tTu-i x [n — 1] 

by 

(ti . . . Tn-2, = T{q, n)Ti . . . Tn-2 = Ti . . . Tn-2T{n, Q ■ Tq-i) 
(n . ..Tn-l)l = Ti . ..Tn-2 

Note that (7r~^)| = (TTj)"^. Note also that we could also define ^ and directly without 
reference to the triangular representation: 

(tt, q)^ = T(n(7r) -|- 1, q)TT* = 7r*T(n(7r) + l,q ■ tt)) 
n = (i"(^^(7i"),g(vr))7r)* = (7rT(/j(7r),Q(7r)))*. 

Proposition 1.1. For n > 2, the maps ^ : 7r„_i x [n — 1] ^ tt^, and {i,q) '■ tTu 
TTn-i X [n — 1] defined above are mutually inverse bijections. Furthermore each restricts to 
C and the restrictions are also mutually inverse bijections. 

Proof. Note that if tt' denotes r(n(7r) -|- 1, q)Tr*, then ^(vr') = q and n{Tr') = n(7r) -I- 1. Thus 
the composition of the maps in either order is the identity. Suppose that a G C„_i and 
1 < q < n — 1. Then p := has the property that n • = • a* for 1 < i < n — 1, and 
hence never equals n. Thus p is an n-cycle. □ 
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By iteration this yields a map iso^ : S„ ^ [n]! := [n] x [n — 1] x ■ ■ ■ x [2] x [1] taking 
a permutation to the sequence of positions j made in the execution of the Fisher- Yates 
algorithm. Sattolo's algorithm fits nicely into this picture. For each i, there is a natural 
inclusion map [i] — > [i + 1]. The product of these gives a map Cn-i : [n — 1]! — > [n]!. 

We summarize the above result in a proposition. 

Proposition 1.2. The following sets are in bijection via the correspondences described 
above. 

(1) the set of possible outputs of the Fisher-Yates algorithm; 

(2) S„; 

(3) the set of triangular decompositions of length n — 1 m S„; 

(4) the set [n] x [n - 1] x • • • x [2] x [1]. 

Furthermore, the measure induced on Sji by the Fisher- Yates algorithm is uniform. 
The following sets are in bijection via the correspondences described above. 

(1) the set of possible outputs of Sattolo's algorithm; 

(2) Cn; 

(3) the set of strict triangular decompositions of length n — 1 in S„; 

(4) the sei [n - 1] X [n - 2] x • • • X [2] x [1]. 

Furthermore, the measure induced on C„ by Sattolo's algorithm is uniform. 

Note that the initial subproduct of length i of a (strict) triangular product ti . . . r„_i 
of transpositions in S„ is itself a (strict) triangular product of transpositions in Sj+i, 
and hence an (i + l)-permutation/cycle. At each stage, forming the next partial product 
involves inserting i + 1 into the current permutation/cycle. This gives an algorithm for 
forming a uniformly random permutation/cycle of a random length; simply form such 
a (strict) triangular product with length chosen according to the desired distribution. 
Clearly the distribution conditioned on the length is uniform. 

Exhaustive generation. There are obvious deterministic versions of the above algo- 
rithm. Instead of randomly choosing the transpositions, we simply run through all such 
transpositions systematically. Every method of generating all elements of [n]l (and the 
corresponding unranking function) can be transferred via the encoding above to a method 
for generating all elements of or Cn+i- 

A common way of enumerating a combinatorial class is to use an incremental method, 
where each object is generated from the last using a small change. The standard minimal- 
change algorithms for permutation generation are given in [Knut2004l 7.2.1.2]. A very 
general method of enumerating permutations is as follows. A Sims table for a subgroup 
G of Sn is a family of subsets Si, . . . of G having the following property: for each j,k 
with 1 < j < k < n, Sk contains exactly one element that fixes all elements greater 
than k and takes k to j, whenever G itself contains such a permutation. It is easily seen 
|Knut20()4l Lemma S] that if 5*1, ... , Sn-i is a Sims table then every element of G has a 
unique representation as a product vr = tti . . .vr„-i, where tt^ G Sk for each k. There is 
also a unique dual representation of the form vr = tt~\ . . . vrf^ with vr^ € 5^, obtained by 
inverting the first representation for vr^^. 

An inspection of the proofs shows that G need not be a group for such results to hold. 
In fact, it is only necessary that G be closed under taking inverses. Thus, for example, 
the set Cn could be used. 

The triangular decomposition fits into this framework. For each k, the set S^ consists 
of all transpositions Tj^ with j < k (in the case S„) or j < k (in the case C„). The Sims 
representation with respect to these sets S^ is precisely the triangular representation. 

For comparison we include the usual inversion encoding in Table [3l 
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Lex order on 


[n] X • 


..x[l] 


Induced order on 




0000 


0100 


0200 


0300 


1230 


3201 


1302 


1203 


0001 


0101 


0201 


0301 


2130 


2301 


3102 


2103 


0010 


0110 


0210 


0310 


2310 


2031 


3012 


2013 


0011 


0111 


0211 


0311 


3210 


0231 


0312 


0213 


0020 


0120 


0220 


0320 


1320 


3021 


1032 


1023 


0021 


0121 


0221 


0321 


3120 


0321 


0132 


0123 



Table 1: The Fisher- Yates encoding for n = 4. 



Lex order on [1] x • • • x 


n 




Induced order on S„ 


0000 0010 0020 0100 0110 0120 

0001 0011 0021 0101 0111 0121 

0002 0012 0022 0102 0112 0122 

0003 0013 0023 0103 0113 0123 


1230 2310 1320 2130 3210 3120 
3201 2031 3021 2301 0231 0321 
1302 3012 1032 3102 0312 0132 
1203 2013 1023 2103 0213 0123 



Table 2: The Fisher- Yates encoding for n = 4. 



Lex order on inversion function 


Induced order on 


0000 0010 0020 0100 0110 0120 

0001 0011 0021 0101 0111 0121 

0002 0012 0022 0102 0112 0122 

0003 0013 0023 0103 0113 0123 


0123 0213 2013 1023 1203 2103 
0132 0231 2031 1032 1230 2130 
0312 0321 2301 1302 1320 2310 
3012 3021 3201 3102 3120 3210 



Table 3: The inversion encoding for n = 4. 



A Gray code for G = S„ is a Hamiltonian path in the Cayley graph of G where the 
generating set is the set of all transpositions. The usual Gray code on words in [1] x • • • x [n] 
induces a Gray code on G via the inversion encoding, since each minimal change to a word 
corresponds to a transposition of adjacent symbols. This can be seen by reading the 
columns alternately downwards and upwards from left to right in Table [3j 

What happens when we use instead the Fisher- Yates encoding? The Gray code order on 
words induces a Hamiltonian path in the Cayley graph of G, but with respect to a different 
set of generators. The generators in question are in fact transpositions and 3-cycles. To 
see this, note that to get from one entry to the next we move from 7rir-7r2 to ttit'tt2, which 
is achieved by multiplying by 7r^^TV7r2. Since r and r' transpose some symbol k with 
respectively, where j ^ j' and < k, the product r'r is the permutation that moves j' 
to j, j to k and k to j' . If either j = k ox j' = k then r'r is a transposition, and otherwise 
it is a 3-cycle. The conjugation by 7r2 preserves the cycle structure. 

Note that the restriction to C is better behaved and the Cayley graph of the set C with 
respect to the set of 3-cycles has a Hamiltonian cycle. Since we always have j, f < A; in the 
case of C, the transpositions are never needed and we always move from one element to 
the next by multiplying by a 3-cycle. Since smaller changes could only be transpositions, 
and multiplying an n-cycle by a transposition can never yield an n-cycle, the enumeration 
described above deserves the name "Gray code for cycles" . 



RANDOM AND EXHAUSTIVE GENERATION OF PERMUTATIONS AND CYCLES 



5 



2. Analysis of some quantities 



Obvious quantities to be studied are: the number of swaps; the number of times a given 
symbol is chosen by the random cahs (we caU this the number of moves, although some of 
these moves will be trivial); the total distance moved by a given symbol; the total distance 
moved. The second and third of these were discussed in |Prod2002l lMahm2003[ IWils2004j 
for the case of Sattolo's algorithm. 

The number of swaps is always n — 1 for each algorithm, but some of these can be 
trivial (and hence executed more quickly) for the Fisher- Yates algorithm, whereas every 
exchange is nontrivial in Sattolo's algorithm. The number of nontrival swaps is the number 
of elements moved by the permutation, or n — /, where / is the number of fixed points. 
The generating function for permutations by size and fixed points is well known to be 



1-Z ' 



For example, the expected number of fixed points is 1 for every n. 

Number of moves and distance moved by an element. To avoid excessive case 
distinctions we consider the slight variant of these algorithms in which the final "swap" 
of 7r[0] with itself is performed (this corresponds to the "downto" loops in Figure [1] going 
down to instead of 1. 

We consider normalized counting generating functions of the form 



^?i(7r) 



7reS,pe[n(7r)] ' '"^"^ ' n>l 

An auxiliary "diagonal" GF will also be useful: 



u 



l<p<n TT£Sn 



G{u,x) := 

TreS 



X 



,n(7r) 



^ n! ^ ^ 

TrSSn 



Here x is a given parameter of interest such as number of moves, etc. Of course F 
and G can be interpreted probabilistically as "grand" PGFs. For example, if xi'^^p) is 
the number of moves made by p in obtaining vr via the Fisher- Yates algorithm, and Mnp 
the random variable obtained by evaluating x element of C„ chosen uniformly at 

random, then letting cj)np{u) = '}Zi>Q^{^np = 0^' denote the PGF of M„p, we have 



F{u,t,x) = Y,X^Y.^''(t>np[u). 

n>l p=l 

We first consider the case where x is the number of moves of a given symbol. The 
triangular decomposition yields the recurrence 



(2.1 



(2.2 



'x{t^I,p) if P 7^ '^(Tr),^ 7^ g(vr); 

1 + X(7ri'9(^)) if P = '^(7r),P 7^ g(vr); 

1 if p 7^ n(7r),p = g(7r); 

^1 ifp = n{7r),p = q{7r). 

In the case where x is the distance moved by an element, we have the recurrence 

'x{iri,p) if p / n(7r),p / g(7r); 

n(7r) - g(7r) + xi-^hqiT^)) if P = n{Tr),p / ^(vr); 
n(7r) — (/(vr) if p 7^ n(7r),p = ^(vr); 

if p = n(7r),p = g(7r). 



x{t^,p) 
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We partition the index set X = {{7r,p) | vr G S, 1 < p < n(7r)} into 4 disjoint subsets 
Xi,...,24 according to the cases just hsted. Denote by 'Ek{u,t,x) the part of the sum 
defining F corresponding to index set Ik, so that F = Si + S2 + S3 + S4. 

Note that for each % we have 



n>l TreSn 

= G{u, tx). 



In the sum Si, indices p € [?^(7r)] satisfying the conditions p 7^ n{7r),p g(7r) occur if 
and only if n(7r) > 2. The set Ii is in bijection with the set 

{(tt, q,p) I n(7r) > 2, 1 < p < n(7r), l<q< n{'!r),p ^ q}. 

Let A{u,t,x) be the antiderivative of F{u,t,x) with respect to x having A{u,t,0) = 0. 
Then for each x we obtain 



Si(.,t,x)=^^^ E ^'-"^'^ 

n>2 7r6S„ l<p<n,p9^q(7r) 



1 



= E^ E E E 

n>2 TreS„-i l<q<n l<p<n,pjtq 

=-EdiijT E *'E«*"" E 

n>l ^ '' l<p<n 7reS„ l<q<n+l,qytp 

= -E7^ E t'E"^""" 

n>l ^ ^ ^' l<p<n TreSn 

= -E^ E " E -Et^ E ('E"*-'" 

n>l ■ l<P<n 7reS„ n>l ^ l<p<n 7reS„ 

= xF{u, t, x) — A{u, t, x). 

We now determine Ti3{u,t,x). The set X3 is in bijection with {tt G S | n(7r) > 2,q ^ 
n(7r)}. Thus for the number of moves we have 

e3(«,<,.) = eSE*""»' = »ES E E*' 

n>2 TreSn n>2 7r|6S„-i q=l 

= > — ^ - = log 1 -ix -ilog 1 -x) . 

n>l 

Note that when t = 1 we have the formula S3(m, 1, x) = «log(l — x) + ux/ (1 — x). 
To obtain S3 for the distance moved, a similar calculation yields 

u 

ulog(l — tx) — tlog(l — ux)] . 



u — t 
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We now consider G. We have for the number of moves 



o«'.-) = E^E»'""' = E^ E + E^ E 

;2 ■ ttgs 

^n(7r|) + l 



n>l ' TTGSn n>2 " 7rGS„,(?(7r)7^n(7r) n>l ' 7reS„,g(7r)=n(7r) 



^ nfyri + 1)! ^ n! ^ 

(7ri,l<9<n(7ri)) 7rGS„,g(7r)=n(7r) 

^n(7r)+l n 

=" E »"""w^)Ti^+«E;j("-«! 

(7r,l<5<n(7r)) \ \ ) ' I „>i 

= nyl('u, 1, x) — « log(l — x). 

Similarly in the case of distance moved we obtain G(u, x) = A{u^ u~^,ux) — log(l — x). 
Thus for the number of moves, by differentiating we obtain the system 

(1 — x)F'{u, t, x) = tG'{u, tx) + T!^{u, t, x) 



G'{u,x) = uF{u,l,x) + 



u 



1-x 

Substituting t = 1 and eliminating G'{u, x) we obtain 

(1 - x)F'{u, 1, x) - uF{u, l,x) = ^— - + S^(u, 1, x) = 
F{u,l,0) = 

which yields 



F{u,l,x) = -—[{l-x)-'-{l-xr^] 
Prom this G can be found explicitly via a single integration. 



G{u, x] ^ 



2 



2-u 



1 , 



+ 



1 — X 1 — u 



— ■ulog(l — x). 



\-u 



To find F explicitly is more difficult, because it requires the integration of (1 — tx) "(1 — 
x)~^, and we do not pursue it here. In any case we have the defining equation 

(1 - x)F'(u, t, x) = utF(u, 1, tx) H h T,'Ju, t, x) 

1 — tx 

F{u,t,0) = 0. 

Using this defining equation we may easily extract the coefficient of x'^t^ to obtain 
the probability generating function 4>npiu), or extract moments by evaluating appropriate 
partial ti-derivatives at u = 1 as usual. For example, (pnpiu) is the coefficient of a;"t^ in 
F, hence equals (l/n) times the coefficient of x'^t^ in xF', and the mean of the random 
variable with PGF (^npiu) can therefore be obtained by evaluating xd'^F/dxdu at u = 1, 
then dividing by n. The defining equation allows us to express these derivatives in terms 
of derivatives of the known series G{u, x) and F(u, 1, x). 

In detail, we see that (with subscripts denoting partial derivatives, and S = E3 to avoid 
notational overload) 

x(l - x)Fis{l,t,x) = txGi2{l,tx) + xT,is{l,t,x) 

Ix 2l-.r- Ixlog(l-lx) Ix 
+ + 7^. - + 



(l-te)2 (l-te)2 1-tx {l-x){l-tx)' 
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Extracting the coefficient of from the right side yields px^ + 2(p — + x^"*"^/ (1 — x) — 
Hp^i where Hp denotes the pth harmonic number X]i<i<p 1/^- Dividing by (1 — x) and 
extracting the coefficient of x" yields (where M denotes the number of moves) 



E[M, 



n + 2p-2- Hp-i 



np\ 



n 



Higher moments can also be obtained with more calculation of the same type, but we 
do not pursue this aspect here. 

We can also immediately extract recurrences for the probability generating function. 
We obtain 

p 

(l>np{u) = -<t>pp{u) + (1 -p/n)u 



u 



n 



1 + XI ^n-l,p{ 
l<p<n-l 



We can now easily write down an explicit formula for the probability generating function 
(pnp by extracting of coefficients from G and using the first recurrence above. We have 

■2 r u{u + l) ...{u + n-2y 



1 



u u 

(Pnn(u) = - + 

n 2 — u ^ 

n-{p-l) p u'^ 

9np{u) = u + 

n n2 — u 

Similarly for the distance moved we obtain 



n! 



u{u + l)...{u + p~2) 



(1 — x)F'(u, t, x) = tG'{u, tx) + T,'^{u, t, x) 

G'{u, x) = uF{u, ,ux) + (1 — x)~^ 

which leads via the substitution t ^ u^^ to 

{1 - x)F'{u,u~^,x) - F{u,u~^,x) = (1 -u"^x)"^ + E^(u,u-\x) 

This equation is exact and leads to 
F(u, u"^ ,ux) = T— ulogfl — x) + Ti^iu, u~^,ux)] 

— -ulogfl — x) + ^ „ [logfl — u^x) — log(l — x)l . 
1 — 

Another integration yields G{u,x) and again we have a defining equation for F{u,t,x). 
Each of these integrations takes us outside the realm of elementary functions. However, 
we can systematically extract coefficients as before from the equations that we have. We 
omit any details of the calculations. The probability generating function ^np of the random 
variable -D„p satisfies the recurrence 

^ n-p 

^np{u) = -Cpp{u) + - 1] 

l + ^U^-P^n-l,p{^) . 
p<n 



1 

1 — MX 



n 
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They have exphcit formulae such as 

1 

n 

Similarly we could extract the mean and higher moments as before. From the explicit forms 
one could consider limit distributions by considering the pointwise limit of the probability 
generating function (pnp or ^„p. We do not pursue this further here as the computations 
are routine but tedious. 

The total distance moved. We consider the total distance moved rightwards by ele- 
ments, Dn = Dnp (this of course equals the total distance moved leftwards by ele- 
ments). As a random variable, -D„ is the sum of -Dn-i and a random variable Ui that is 
uniform on [0..n — 1]. Thus Dn is distributed as Yl^=2 probability generating 

function nr=iT^- 

Note that this PGF is the same as the one for inversions. Thus the number of inversions 
and the total rightward distance have the same distribution (in otlier words, Dn a- 
Mahonian statistic). Hence for each A:, the number of permutations in S-n, with h 
inversions is the same as the number of permutations in whose rightward distance is 



n-l 



in-3 



+ 



l-V? 



3. Extensions and discussion 

Despite an extensive literature search, I can only find two places in the literature in 
which the very natural Fisher- Yates encoding is mentioned. In neither paper was the 
connection with cyclic permutations mentioned. 

What we have called the Fisher- Yates encoding was used in |MaRa2001| to study anti- 
excedances. In |MyRu2001| , two unranking and ranking functions for permutations were 
presented, each taking linear time to compute. Although not mentioned in that paper, it is 
easily seen that those orderings correspond via the Fisher- Yates encoding to lexicographic 
order on the triangular cartesian product [1] x • • • x [n] or [n] x • • • x [1] as we have described 
above. 

We note that Sattolo's algorithm is a special case of a method to uniform generation of 
permutations with a fixed number k of cycles |Wilf| . 

Suppose that m > is fixed, and at the ith step, we swap 7r(z) and 7r(j) where j is 
chosen uniformly from [l..n + 1 — i — m\. The cases m = and m = 1 respectively 
correspond to the Fisher- Yates and Sattolo algorithms. Other values of m do not appear 
to be particularly interesting, although we have not pursued them. 

I have not yet been able to determine whether the statistic D above is known. An 
interesting question is its correlation with other well-known permutation statistics. 

There is a small connection between the Fisher- Yates algorithm and sorting. 

Proposition 3.1. Let vr E S„ and let ti . . .r„_i be its triangular decomposition. Then in 
order to sort Tr~^ , selection sort applies the transpositions Tn-i, • • • ,ti in that order. 

Proof. Selection sort first chooses n and puts it in the correct position; this corresponds 
to postmultiplication by T{n,mT), which corresponds to premultiplication by r(n, nvr"^). 
The result follows by induction. □ 

For example, to generate vr = 3421 the algorithm proceeds as follows: 1234, 4231, 4321, 
3421, yielding the strict triangular decomposition (12)(23)(14). Applying these in turn 
to TT~^ gives 4312, 2314, 2134, 1234 which is the list created by selection sort when sorting 
vr"^ = 4312. 
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